BI-TP 2006/05 
RBRC 597 

Transverse flow in relativistic viscous hydrodynamics 

Rudolf Baier/ Paul Romatschke/ and Urs Achim Wiedemann^' ^ 

^ Fakultdt fiir Physik, Universitdt Bielefeld, D-33501 Bielefeld, Germany 

^Department of Physics and Astronomy, University of Stony Brook, NY 11794, USA 

^RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 

(Dated: February 9, 2008) 

Hydrodynamic model simulations of Au-Au collisions at RHIC have indicated recently, 
that with improved simulations in the coming years, it may be feasible to quantify the vis- 
cosity of the matter produced in heavy ion collisions. To this end, a consistent fluid dynamic 
description of viscous effects is clearly needed. In this note, we observe that a recently used, 
approximate form of the 2nd order Israel-Stewart viscous hydrodynamic equations of motion 
f^ \ cannot account consistently for the transverse flow fields, which are expected to develop in 

^SJ ' heavy ion collisions. We identify the appropriate equations of motion. 
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I. INTRODUCTION 

Ultra-relativistic nucleus-nucleus collisions exhibit large collective phenomena, such as trans- 
verse flow. In the discussion of the dynamical origin of this collectivity, relativistic hydrodynamics 
plays a central role. This is so, since the densities attained in nucleus-nucleus collisions imply very 
^^ , small mean free paths (< 1 fm), a condition, which hydrodynamics can naturally account for, while 

Cp I other approaches to multi-particle dynamics (such as parton cascades) have difficulties to control. 

(<— ^ ■ Indeed, simulations based on ideal fluid dynamics successfully describe the main characteristics of 

\C) • soft momentum distributions in nucleus-nucleus collisions at RHIC [l|, y,y{. However, the results 

^^ \ of these simulations may depend significantly on the modeling of the initial conditions and freeze- 

out. This raises the question whether the success of ideal hydrodynamics is indicative of an almost 
viscosity- free matter produced at RHIC, or to what extent matter of non- negligible viscosity can 
O \ be accounted for by ideal hydrodynamics due to a compensatory choice of initial conditions and 

freeze-out jj]. To address this point, the study of the hydro dynamical equations of motion in the 
presence of viscous corrections is clearly needed. 

Including dissipative corrections of ideal fluid dynamics to first order in a gradient expansion is 
K> \ known to lead to an acausal dynamics, exhibiting significant unphysical features |5|. A consistent 

j^ ' relativistic dissipative hydrodynamics requires a gradient expansion to second order [a, lD| • The 

application of this Israel-Stewart theory to relativistic heavy ion collisions is at the very begin- 
ning |J, la ISl; IliJ) lUJ, U2] • III this note, we show that a specific approximation of the Israel-Stewart 
equations of motion [ll|, [l^] , used recently, does not allow for a consistent description of transverse 
flow, and we identify the full equations of motion, which should be used. 

For simplicity, we consider hot matter with an ultra-relativistic equation of state (EoS) e = 3p. 
In the absence of conserved charges (vanishing chemical potential) and for finite shear viscosity r/, 
the energy-momentum tensor takes the form 

T^^ = (e + p)u^'v!' - pg^"" + U^" . (1) 

Here €,p are the energy density and pressure, u'^ denotes the local fluid velocity, u^u^ = 1, and 
H^'^ is the traceless shear tensor. The tensor T^'^ is conserved, 

d^T^'''=0, (2) 

where d'^ is the covariant derivative. In the Israel-Stewart theory, the evolution equation of the 
shear tensor is given by |j, |3| 

TnA^A^DH"^ + W^ = rjiV^'u") (3) 



— n^.n^'^>o. (4) 



with V = A^'^d^, A'^'^ = 5^^^ - u^'u", {V ^u^) = V ^Uy + V^u^ - |A^^V„u", and D = u^'d^. 
In the limit of vanishing relaxation time rn, one finds the defining equation for 11'^'' in first order 
dissipative hydrodynamics. Here and in what follows, we neglect for simplicity other dissipative 
corrections (bulk viscosity, heat conductivity), as well as contributions from the vorticity to ©• 

Eq. © can be derived from kinetic theory, which is one way of determining all possible terms 
entering the equation of motion jj|. An alternative 'short-cut' derivation is based on the constraint 
that up to second order in a gradient expansion, the comoving entropy density s^ cannot decrease 
0,0 

Td^s^ = n^.i-rnl^n^" + r? < V^'u'' >]/(27?) 
1 

2ri' 
Motivated by this constraint, one uses in the recent literature jlll. Il2| the equation of motion 

mDIi^'' + n'^'^ = r? < V^'u'' > . (5) 

However, since n^H'*^ = 0, Eq. Q does not specify the components parallel to u^ in ©. Thus, 
Eq. ((21) is based on @ and the additional approximation that u^DU'^'^ = 0, which one obtains 
by contracting Q with n^. On the other hand, contracting the full Israel-Stewart equation Q 
with Ufj_, one finds u^DIi^'^ = —U'^'^Du^, which is not an additional constraint but a simple 
consequence of the orthogonality u^Ii'^'^ = 0. It follows that involves the additional assumption 
U^'^DUfi = 0. At least for the physically relevant case of a radially symmetric collision region, this 
latter constraint implies necessarily (see the appendix for details of this argument) 

Du'' = 0. (6) 

We emphasize that © is a necessary consistency condition for the approximate equation of motion 
©, but it is in general not satisfied for the full Israel-Stewart equation of motion Q. 

II. SOLUTION FOR IDEAL HYDRODYNAMIC EVOLUTION 

Is the additional assumption Q acceptable for the simulation of relativistic heavy ion collisions? 
To address this question, we study here the resultin g so lution of ideal hydrodynamics, which serves 
as the reference in comparison with viscous effects [ul, |l3; llJ, Il5| • Our starting point is Eq. ©, 
which in the sense of a 'proof by contradiction' is assumed to be valid in this section. We also use 

^^^ = 0, (7) 

De = -ie + p)V^u^' ^ Dp = --p V^u^^ , (8) 

which follow from the conservation laws ©, together with ©• We consider a system with longitu- 
dinal boost-invariance [1^] and cylindrical symmetry in the transverse direction. The independent 
variables are r, r, cp and rapidity t]. Because of the symmetries, the pressure p = p{t, r) and radial 
fiow velocity f = u(r, r). 
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v = — ,v<l, (9) 

are functions of r and r only. Using the metric components, grr = l;5rr = —^,9rjr] = —T'^,g<i><j) = 
— r'^, the relativistic hydrodynamic equations become 

{dr + vdr)p{T,r) = 0, (10) 

{dr + vdr)p{T,r) + -p{T,r)[-f'^{dr+vdr)v{T,r) + - + -] =0, (11) 
6 TV 

{dr+vdr)v{T,r) = 0, (12) 



dr = ^- This set of quasilinear partial differential equations 



where 7"^ = 1 — v^, and dr — -q;^, 

can be solved by the method of characteristics. 

Prom Eq. (|12() . one finds that along the characteristic T,r{T) that fulfills the equation -^ =v, 
the value of v is constant, 



dv 



0. 



Therefore, one can integrate J dr = J vdr, obtaining 

VT — r = const = f{v) 



(13) 



(14) 



where we used jr = to rewrite the constant as an arbitrary function of v. Eq. ()14|) is the general 



dr 
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solution of Eq. (fT2|) . It has been discussed previously in a different context, see e.g. Biro 17l.ll8||. 
We now use 



^'^{dr + vdr)v = drV, p = p (vt) 



-4/3 



to rewrite Eq. (|11() as 



{dr + vdr)p = --pdrV. 
o 

Using again the method of characteristics, we find 

d Inp ' = —{dr-v) dr, vdr = dr. 



(15) 



(16) 



(17) 



Since the second equation is the same characteristic as before we know that along this characteristic 
^ = 0. Rewriting 

1 



dr V 



r - f'{v) 



from Eq. 1)14(1 with f'{v) = -J-, one can thus integrate Eq. (|17() . obtaining 



p-3/4 = (^ _ /'(„)) X const. = (r - nv))g{v), 



P 



(18) 



(19) 



where we introduced another arbitrary function g{v) for the same reasons as above. In summary, 
we find for the general solution of Eq. (|lll) 



P 



m = ^rTr\drv)-g{v) , 



(20) 



where g{v) = g ^{v). The remaining Eq. (|10|) can be used to constrain the functions f{v) and 
g{v). Inserting Eq. ^1 into Eq. (UHl) and using Eqs. (|12I14I18|) we find 

-vfgif'f 

+r [(1 - v^)f f'g' + g{l + v^){ff + gfiSvf + (1 - v^)/") 
(1 - v^)if + vf')g' + g (2vf + (2 + St;^)/' + v{l - v^)f") 



+ T^ 



il + 2v')g + vil-v')g' 



0. 



(21) 



For solutions relevant in heavy ion collisions, we require that at all spatial points r and for all 
times T > 0, the pressure and its time derivative are finite. In particular, we require < p(r, r = 
0) < 00, and I drp{T,r = 0) | < 00. According to Eq. (|2U() . the first condition leads either to (i) 



liuir^idrV ~ r or to (ii) lim.r<^ig{v) ~ r. In the next paragraph, we shah show that case (i) 
is reahzed. Because of Eq. (|12|) . case (i) impHes drv{T,r = 0) = 0. For case (ii), one considers 
drv\r=o / 0. Then, to satisfy drP < oo, one has to require hmr<i drg{v) ~ r, and this leads again 
to drv{T, r = 0) = 0. As a consequence of Eq. (fTl|) we thus find in both cases v{t, r = 0) = 0. 

Since v{t, r = 0) = 0, obviously v does not depend on r at r = 0. Thus, for Eq. (PT|) to hold 
at all times r, each power of r in (|2H) has to vanish separately. One can exploit this property in 
a neighborhood around ?; = by doing a Taylor expansion of (|21j) in v around v = 0. One then 
finds that not only the u-dependent prefactors of the powers r" in (|21() must vanish at u = 0, but 
also all their derivatives must vanish. Thus, these prefactors must be zero and one finds 

(22) 





f{v) = 0, 


g{v) = const x -(1 — v"^)^'"^ 


which leads to 








r 
v{T,r) = -, 

T 


p{T,r) = const X (r — r )' 



(23) 

This velocity field v and pressure p are physical for r < t, only. At r = the pressure behaves as 
p{T,r = 0) ~ I/t'^, which differs significantly from the Bjorken cooling law p{t) ~ t~'^'^ UM- This 
exact analytic solution (|2,3|) of ideal relativistic hydrodynamics in (1 + 3) dimensions, constrained 
by Du^' = and restricted by boost invariance and cylindrical symmetry, has also been discussed 
previously in jlfll |. 

For hydrodynamic simulations of central heavy ion collisions, the initial conditions satisfy usu- 
ally v{To,r) = for finite time tq and all values of r UJ, Il2l . Il5l |. This initial condition is in- 



compatible with the requirement ()22() that f{v) = in a finite region around v = 0, and so it 
is incompatible with Eq. (|21() . Moreover, in numerical solutions of the full ideal hydrodynamic 
equations of motion (for which Du^ 7^ in general), one finds that the velocity field for small 
but fixed values of r increases with r for some time (see e.g. Fig. 1 in |l5|). In contrast, for 
the solution (jSSI), the velocity v{T,r = const) decreases with r. This illustrates, that Du'^ = 
is not a physically justified approximation in the context of heavy ion collisions, and so it should 
not be adopted for the refined formulation of viscous hydrodynamics. For reliable results, the full 
Israel-Stewart equation © rather than Eq. © has to be used. 
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APPENDIX A: EQUATION © IMPLIES Du^' = 

In this appendix, we consider a longitudinally boost-invariant system with cylindrical symmetry 
in the transverse plane. We show that approximating the full hydrodynamic equations © by © 



amounts to the assumption Du^^ = 0. We start from the observation that Eq. Q agrees with 
Eq. ®, if 

WDua = 0. (Al) 

Since only u^,u^ are non- vanishing, we find (choosing ^ = r) 

W^Du'' - WDvr = 0. (A2) 



Using UylV^^ = and vT' = yjl + (u**)^, this expression can be written as 

2 \ 




- 1 LDu'^ = 0. (A3) 

Since u^ is in general unequal u"^, this means that Du^^ = if H^^ is non- vanishing. In contrast, 
for the full viscous equation of motion 0, Du^^ 7^ in general. 
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